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Differential cross sections for the reaction j p ^ p tv^'' have been measured with the CEBAF 
Large Acceptance Spectrometer (CLAS) and a tagged photon beam with energies from 0.675 to 
2.875 GeV. The results reported here possess greater accuracy in the absolute normalization than 
previous measurements. They disagree with recent CB-ELSA measurements for the process at 
forward scattering angles. Agreement with the SAID and MAID fits is found below 1 GeV. The 
present set of cross sections has been incorporated into the SAID database, and exploratory fits 
have been extended to 3 GeV. Resonance couplings have been extracted and compared to previous 
determinations. 
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I. INTRODUCTION 

The spectrum of baryon resonances has been exten- 
sively explored for clues to the internal structure of nu- 
cleons. Experimental and phenomenological programs 
are working in tandem to refine and expand the known 
resonance masses, widths, and electromagnetic couplings, 
which provide tight constraints for QCD-inspired models, 
and valuable benchmarks for lattice calculations. The 
most precisely determined resonance properties are asso- 
ciated with low-lying states and those higher-mass states 
with clear Breit-Wigner signatures in amplitudes where 
they contribute. Many other states have been observed, 
but with widely- varying mass and width estimates. 

Much of the non-strange baryon spectrum has been de- 
duced from fits to pion-nucleon scattering and photopro- 
duction data. However there are difficulties encountered 
in extracting the states from experimental data. While 
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many states have been predicted by QCD-inspired mod- 
els, far fewer have been clearly identified. Some of these 
missing states may be weakly coupled to the pion-nucleon 
channel, and this possibility has motivated the study of 
other channels (involving 77 TV, KA, and ttttN final states, 
for example). Many of the states established in fits to 
elastic pion-nucleon scattering data have ttN branching 
ratios jj only of order 10-20% and therefore are not ex- 
pected to be easily seen in irp elastic scattering. This is 
particularly true for resonances with masses above about 
1.7 GeV. 

Photodecay amplitudes can be extracted from fits to 
the available database of charged and neutral pion pho- 
toproduction data. A knowledge of the resonances con- 
tributing to pion-nucleon elastic scattering is crucial in 
this task. Precise measurements of the cross section and 
polarization asymmetries are equally important. How- 
ever, until very recently, such measurements were lack- 
ing in the region above 1 GeV in photon energy. Recent 
Jefferson Lab measurements have shown that stirprises 
are possible in the behavior of both polarized [2| and 
unpolarized 3 cross sections. Present fits covering pho- 
ton energies to 2 GeV and beyond are necessarily model- 
dependent and undcrconstraincd. 

In this paper, we report measurements of the unpolar- 
ized differential cross sections of neutral pion photopro- 
duction on the proton for incident photon energies from 
0.675 to 2.875 GeV. We have included the present set 
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of CLAS cross sections in a multipole fit to the avail- 
able data covering the resonance region. At the highest 
energies, recent CB-ELSA Q measurements are avail- 
able, and we have compared both data and fits to show 
where deviations occur. Resonance couplings have been 
extracted for those states that give a significant contri- 
bution to the photoproduction process, fn some cases, 
these couplings differ significantly from previous deter- 
minations. 

The paper is laid out in the following manner: We give 
a brief overview of the experiment in Section [TTl A more 
detailed examination of methods used in the data reduc- 
tion follows in Section [ml The uncertainty estimates for 
the cross sections obtained are given in Section IIVI The 
experimental results are described in Section |Vl Various 
fits to the data are described in Sectioi TvH and the under- 
lying multipole amplitudes and resonance contributions 
are displayed and compared to previous determinations 
in Section lVIII Finally, in Section lVIIIi we provide a brief 
summary of the results of this study and consider what 
extensions of this work would be particularly helpful in 
the future. 



II. EXPERIMENT 

The differential cross sections for the reaction 7 p — > 
p tt" were measured with the CEBAF Large Acceptance 
Spectrometer (CLAS) Q and the bremsstrahlung pho- 
ton tagging facility ( "photon tagger" ) @ in Hall B of the 
Thomas Jefferson National Accelerator Facility (JLab) 
as part of a set of experiments running at the same time 
with the same experimental configuration (cryogenic tar- 
get, tagger, and CLAS) called the "glc" run period. The 
cross sections can be found in electronic form in Ref. Q 
and were part of a program of meson photoproduction 
measurements undertaken using the CLAS and photon 
tagger g, d, [13, Hi!, d, [H, Q [H . The database entries 
include the differential cross sections, as well as uncer- 
tainties (excluding the overall absolute normalization un- 
certainty), for each incident photon energy and cos 
shown in this document. 

The full data described here consisted of two running 
periods with different incident electron beam energies of 
2.445 and 3.115 GcV. Tagged photons were incident on 
an 18-cm-long liquid hydrogen target placed at the center 
of CLAS. This target was enclosed by a scintillator array 
(start counter) that detected the passage of charged par- 
ticles into CLAS from the target ^T§\ . The event trigger 
required the coincidence of a post-bremsstrahlung elec- 
tron passing through the focal plane of the photon tag- 
ger and at least one charged particle detected in CLAS 
and the start counter. Tracking of the charged particles 
through the magnetic field within CLAS by drift cham- 
bers provided determination of their charge, momentum, 
and scattering angle. This information, together with 
the particle velocity measured by the time-of-flight sys- 
tem [l7| and start counter, provided particle identifica- 



tion for each particle detected in CLAS and its corre- 
sponding momentum four-vector. The methods used for 
extracting the differential cross sections for 7r° photopro- 
duction are presented in the next several sections. The 
technique is outlined initially, and then each step is de- 
scribed in further detail, with data and tests that support 
the validity of the approach. 



III. DATA REDUCTION 

Yields for the tt^ mesons resulting from photoproduc- 
tion on the protons within the cryogenic hydrogen target 
were determined using the missing mass technique, as- 
suming the two-body reaction 7 p — > p A", where X is 
the particle hypothesized to be missing. Briefly, the fol- 
lowing steps are involved in this analysis: 

1. Identify the recoil proton in CLAS, determin- 
ing the scattering angle and momentum for 
the proton. 

2. Calculate the missing mass for the recoil pro- 
ton based on the assumption that the reaction 
observed is 7 p ^ p AT. 

3. Distribute the events of the resulting miss- 
ing mass spectra into bins defined by incident 
photon energy and tt*^ center-of-mass scat- 
tering angle cos. 6^^ , based on the assump- 
tion that the reaction observed is 7 p — > p tt*^. 

4. Identify the tt'^ meson missing mass peak. 

5. Determine the yield for the tt'^ meson by sub- 
tracting any background from beneath the 
meson peak. 

6. Correct the meson yield for spectrometer ac- 
ceptance and detection efficiency based on 
Monte Carlo estimates of those quantities. 
This Monte Carlo technique was compared 
to empirical measurements of the CLAS ac- 
ceptance and detection efficiency for a single 
proton, using measurements of the reaction 
7 p — > p TT^ TT'^ made concurrently with the 
data discussed here, to validate the Monte 
Carlo simulation used. 

7. Normalize the yield using a measured abso- 
lute photon flux normalization procedure. 

In the following sections, each of these steps is de- 
scribed. Also presented are sample results, and, in some 
cases, tests that establish the validity of the procedures 
used. 



4 



CO 












o 






P 








: J 








\ K 














V fl j" 










1 D ;- 









Mass !G^>' ' c") 



FIG. 1: Mass distribution of the detected charged hadrons 
(note logarithmic scale). 



A. Recoil proton and pion identification; kinematic 
variables 

The tracking information provided by the drift cham- 
bers within CLAS gave momentum and scattering angle 
information on charged particles scattered within the de- 
tector volume. Time-of-flight and start counter informa- 
tion, coupled with the track information provided by the 
drift chambers, yielded velocity and momentum determi- 
nation. 

Particle identification in this analysis was performed 
using the GPID algorithm [l^ . The method uses the mo- 
mentum of the detected particle, and sequentially calcu- 
lates trial values of /? for the particle for possible particle 
identities. Each one of the possible identities is tested by 
comparing the trial value of P for a given particle type 
to the empirically measured value of (3 (as determined by 
CLAS tracking and time-of-flight information) . The par- 
ticle is assigned the identity that provides the closest trial 
value of (3 to the empirically measured value of (3. Figure 
[T] shows the mass distribution of the identified positively 
charged particles. The GPID algorithm also attempts to 
find a matching photon in the tagging system for every 
charged particle detected in CLAS. A matched photon 
means that there was one and only one tagged photon 
in the trigger window, which, in this analysis, is defined 
as being within the trigger coincidence window. Parti- 
cles that are determined not to have a matching photon 
are considered to be a measure of the accidentals (to be 
described in more detail in the next subsection). 

Geometrical fiducial cuts in each of the six sectors of 



CLAS were imposed on all protons. The region selected 
for accepting protons in each sector corresponded to a re- 
gion of relatively uniform detection efficiency versus az- 
imuthal angle. 



B. Missing mass reconstruction 

The momentum determined by CLAS was corrected 
for energy loss within the cryogenic target cell to recon- 
struct the momentum of the detected proton at the re- 
action vertex within the cryogenic target. In addition 
to the energy loss correction, a further momentum and 
photon energy correction developed by Williams et al. 

for glc data was applied. (This second momentum 
correction was required to correct for problems with the 
magnetic field map associated with the CLAS detector.) 
The measured scattering angle and momentum can be 
used to construct a missing mass based on the assump- 
tion that the reaction observed is 7 p p X, where 
X is the undetected particle in the two-body final state. 
Based on this assumption, the missing mass spectrum of 
the full spectrometer acceptance for all photon energies 
is shown in Fig. [21 The 77, p'^/t^J, and 77' peaks are clearly 
seen atop a background dominated by multi-pion events. 
The 7r° peak is clearly discernible. 

Taking each proton event that did not have a match- 
ing incident photon as noted above, and integrating over 
all of the out-of-time (not within the trigger coincidence 
window) incident photons for that event, determined the 
distribution of accidental coincidences between the CLAS 
and the photon tagger, under the assumption that cou- 
pling the out-of-time tagger events for each non-matched 
proton created a fair representation of the accidental co- 
incidences between the CLAS and tagger. 

A failure to match a particle to an incident photon 
mainly occurs when reconstructed timing information for 
the track or tagging system is missing. 



C. Distribution of events into kinematic bins 

The events, from the 2.445 and 3.115 GeV data sets, 
constituting the full missing mass spectrum described in 
the previous section were distributed into bins in photon 
energy and cos6'^^ . The widths of these "kinematic 
bins" (A_B-y = 50 MeV in photon energy and Acos6'Jj„ 
= 0.1) were chosen such that, in general, there were at 
least 1000 tt" events in each kinematic bin. 



D. Determination of pion yield within each 
kinematic bin 

The 7r° yield within each kinematic bin was determined 
by subtracting the background under the peak in the 
missing mass spectrum. We proceeded with the assump- 



5 




0.6 O.S 1 1.2 1.4 

Missing Mass (GeV/c^) 




-0.04 -0.02 



0.05 0.08 0.' 

MM' (GeVVc*) 



FIG. 2: Missing mass spectra obtained using CLAS, assum- 
ing the reaction 7 p ^ p X. 

tion that the background in the missing mass spectra 
arose from two particular sources: 

1. Accidental coincidences between CLAS and 
the photon tagger. 

2. Two pion photoproduction via the reaction 
7 p — > p X, where X = tt+tt^. 

The accidental contributions to each kinematic bin were 
determined as described in Section III.B. Since events 
with pTT+TT^ final states were copiously produced in this 
experiment (as may be seen in Fig. [2]), the shape of the 
low energy portion of the 27r background contribution 
could be reliably determined. This shape was used to 
generate the background beneath the n'^ peak, which was 
then subtracted from the tt'^ yield for each kinematic bin. 
An example of the individual contributions to the back- 
ground beneath the pion peak can be seen in Fig. [31 

E. Acceptance and efficiency 

The CLAS detector acceptance and detection efficiency 
for recoil protons were measured and then compared to 
a Monte Carlo simulation. The reaction used to deter- 
mine the empirical acceptance and efficiency of protons 
in CLAS was 7 p — > p 7r~ tt"*". This reaction was used 
due to the high number of events for that final state and, 
since all the final products leave charged tracks in the 
CLAS, these events are easily observed. Both pions were 
required to be detected in the event and both matched 
to the same photon. 



FIG. 3: n meson yield extraction for the kinematic bin with 
E-i = 1.425 GeV and cos6lc.m. = 0.45. The top panel a) is 
the missing mass yield for this bin, with the accidental con- 
tribution displayed as a shaded region. The middle panel b) 
shows the missing mass distribution with the accidentals sub- 
tracted, and the shaded region represents the 27r contribution. 
The bottom panel c) shows the extracted tt" yield after both 
contributions have been subtracted from the missing mass 
distribution. 



The data used for the empirical acceptance and effi- 
ciency calculations included only events where two and 
only two charged pions were detected in the CLAS. For 
each event, a missing mass reconstruction from the kine- 
matical information from the two pions was performed 
to determine if a proton should have been seen in the 
CLAS. As shown in Fig.[31 the proton generally was very 
cleanly defined, so that the determination that a proton 
should have been seen could be made without ambiguity. 
The background beneath this peak was approximated by 
a third order polynomial and subtracted. The same fidu- 
cial cuts applied to the protons noted above were applied 
to both reconstructed and CLAS-identified protons. 

The proton acceptance and efficiency a empirically de- 
termined in this fashion is the ratio of the number of 
protons "seen" Ng (i.e., identified in CLAS through the 
normal particle identification procedure, GPID, in the 
empirical data set) and the total number of events Np 
where a proton should have been seen, based on recon- 
structed four-vectors from the pions, detected in the same 
subset of events. 

A comparison of Monte Carlo events to actual data for 
the 7 p — > p 7r~ 7r+ reaction (rebinned for the 7 p — > p Tr*^ 
reaction) was performed. Simulated events were obtained 
by generating 10^ 7 p — > p Tr~ tt"*" events that were 
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FIG. 4; Missing mass (mass of X) for the reaction 7 p — > 
TT^ TT^X near the mass of the proton (shaded region repre- 
sents the background). 



isotropic in phase space and then processed through a 
fuU GEANT simulation of CLAS. These events were then 
"smeared" to simulate the drift chamber and time-of- 
flight resolution. After smearing, the events were pro- 
cessed as normal data. 

A comparison between the empirical and Monte Carlo 
data can be seen for an incident photon energy of 725 
MeV in Fig. [5] To perform a statistical comparison for 
these two data sets, the bottom panel of Fig. [5] (Monte 
Carlo/empirical) was fit with a polynomial of zero order 
over the cos0'^^ range of -0.85 to 0.45. For angles above 
this range, there were not enough empirical events, af- 
ter rebinning to the ir^ reaction, to perform a reliable 
comparison. For cosd'^^ < —0.9, the Monte Carlo did 
not agree well with the empirical method, and, for this 
reason, points with cos0^^ < —0.9 were left out of the 
subsequent steps of analysis. The zero-order result of 
the fit is called the "acceptance and efficiency ratio." 
The ratio of empirical to Monte Carlo for each energy 
and cosOl^ bin (within 0.65 GeV < < 1.8 GeV and 
—0.9 < cosO'^^ < 0.5) was placed in a histogram, and 
then fit with a Gaussian. The center of the Gaussian 
was found to be 1.0012 ± 0.002, with standard deviation 
0.0312 ± 0.016, and reduced — 0.76, demonstrating 
our ability to find and correct for inefficiencies and ac- 
ceptance. 

With this confirmation of the Monte Carlo validity, 
for -0.9 < cose/^.^ < 0.5 and 0.65 GeV < < 1.8 
GeV (and assuming the validity for all other energies and 
—0.9 < cosd^^ < 0.9), the acceptance and efficiency 



FIG. 5: Empirical and Monte Carlo acceptance and efficien- 
cies a for Tv^' photoproduction for Ej = 725 MeV. The top 
panel shows a (filled circles represent the empirical method 
and open circles represent the Monte Carlo method). The 
bottom panel shows the ratio of a (Monte Carlo method di- 
vided by the empirical method). 



for each kinematic bin for the reaction 7 p — > p tt'^ was 
obtained by generating 10'' events (weighted by the cross 
sections given by the SAID solution [205). The ratio of 
processed events, to the number generated, for a given 
bin, served as a measure of the acceptance and efficiency. 



F. Sector-by-sector comparison 

The CLAS has six sectors that are nominally iden- 
tical, but these sectors, in fact, differ relative to each 
other in acceptance and detection efficiency due to dif- 
ferences that either initially existed in their construction 
or by various changes and differences that have arisen 
over time since they were first installed. A sector-by- 
sector comparison of the differential cross sections was 
performed to check the consistency of the extracted cross 
sections. Since the Monte Carlo simulation should reflect 
sector-by-sector changes in the detector arising from, for 
example, broken drift chamber wires and time-of-flight 
paddles, a sector-by-sector comparison of the differential 
cross section explores the reliability of the Monte Carlo 
simulation with respect to these detector irregularities. 
A typical sector comparison plot is shown in Fig. [6l The 
plot shows the differential cross sections of each sector 
of CLAS divided by the average over all sectors. An 
examination of these plots for each energy revealed no 
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tailed discussion of the g2a experiment can be found in 
Ref. [2l|.) By observing events in the g2a running pe- 
riod that had a photon trigger and no charged trigger, 
yet clearly had a proton detected by CLAS in the event, 
the inefficiency of the charged particle trigger in CLAS 
for protons could be determined. This correction was 
applied to each kinematic bin and was always less than 
1.0%. 
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FIG. 6: Differential cross sections for each sector divided by 
the average over all CLAS sectors {E~, — 1.025 GeV) for the 
reaction j p ^ p ■k'^. 



systematic shifts visible in the differential cross section 
for either or cos 61!^ . 



In contrast to many of the previously published data 
sets for 7 p — > p 7r°, the data in this analysis were not 
normalized to previous data for the process or to the 
SAID solution for any reaction channel. The absolute 
photon flux for the entire tagger photon energy range 
was determined by measuring the rate of scattered elec- 
trons detected in each counter of the focal plane of the 
bremsstrahlung photon tagger by sampling focal plane 
hits not in coincidence with CLAS. The detection rate 
for the scattered electrons was integrated over the life- 
time of the experiment and converted to the total number 
of photons on target for each counter of the tagger focal 
plane. The tagging efficiency was measured in dedicated 
runs with a total absorption counter (TAC) downstream 
of the cryogenic target, which directly counted all pho- 
tons in the beam. The details of the method can be found 
in Ref. 



G. Bin migration 

To estimate the systematic error associated with bin 
migration, the acceptance and efficiency results calcu- 
lated using SAlD-weighted events were compared to ac- 
ceptance and efficiency results using non-weighted events. 
For each kinematic bin, the systematic error associated 
with bin migration would be the uncertainty in the cor- 
rected value multiplied by the amount of the correction. 
Since the amount of the correction was found to be typ- 
ically less than 1.0%, and always less than 2%, the sys- 
tematic uncertainty associated with bin migration was 
assumed to be ignorable. 



H. Trigger inefficiency 

The determination of a charged particle trigger inef- 
ficiency for the data was performed by looking at data 
from a running period that had a CLAS-detected photon 
event trigger condition in addition to the CLAS-detected 
charged particle event trigger described above for the glc 
running period. An experimental running period that 
had both charged particle and photon triggers occurred 
just before the running period used in this analysis. (This 
dual trigger running period was called g2a; a more de- 



IV. UNCERTAINTIES 

Having investigated various sources of uncertainties in 
the analysis, we collect and summarize here the various 
uncertainties determined for the cross sections obtained 
in this work. 



A. Trigger inefficiency 

From the trigger inefflciency study given in subsection 
III.H, an overall estimated systematic uncertainty of 1% 
for the trigger inefficiency correction factors was taken 
as a very conservative estimate of the systematic uncer- 
tainty associated with the trigger inefficiency. 



B. Baclcground subtraction 

The uncertainty associated with the background sub- 
traction is purely statistical, and these were taken into 
account on a bin-by-bin basis. No systematic uncertain- 
ties for the background subtraction are included. 
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C. Effects of momentum correction 

As noted in Section III.B, a pair of momentum correc- 
tions for the rescattered proton was made prior to form- 
ing the missing mass spectra. These momentum correc- 
tions affect extraction of the 7r° from the missing mass 
distribution, and can alter the the center-of-mass angle 
for the scattered proton, sometimes transferring events 
from one kinematic bin to another ( "bin migration" ) . An 
estimate of the systematic uncertainty introduced by ef- 
fects of these momentum corrections can be formed by 
looking at the magnitude of the momentum correction 
and the magnitude of the resulting deviation from the 
optimal value determined by a pull distribution. 

To determine the uncertainty associated with the mo- 
mentum correction, the reaction ^ p p ■k~ 7r+ was 
studied using methods described in greater detail in 
Ref. 19]. The variable z is defined as z = /S.pp/a = 
{Pp(fit) - Ppimeasured)) /<^, whcrc Pp(/jt) represents the 
best value of the momentum as determined by a kine- 
matic fit for data which does not include the CLAS mea- 
sured value for the proton momentum, and a is the stan- 
dard deviation. The distribution of z was examined by 
histograming and fitting the results to a Gaussian. His- 
tograms were generated for z from events with and with- 
out the momentum correction for the proton. The differ- 
ence between these means is taken as a measure of the 
average momentum correction. 

The initial (final) mean and standard deviation of the 
pull distribution were 0.208 and 0.963 (-0.022 and 1.011). 
The magnitude of the change in the mean before and af- 
ter the momentum correction is 0.23. The amount that 
the mean of the distribution is still different from the 
optimal value of zero is 0.022. This suggests that the 
momentum correction is only good up to 0.022/0.23 — 
0.096. Thus, we estimate conservatively that the sys- 
tematic uncertainty in the momentum correction is 10%. 

The resulting uncertainty in the cross section due to 
the uncertainty in the momentum correction is taken to 
be 10% of the absolute value of the difference between 
the momentum corrected and non-momentum corrected 
differential cross sections. This uncertainty is typically 
less than a few percent and is added in quadrature on a 
bin-by-bin basis. 

D. Photon flux normalization 

Due to the large number of out-of-time photons used to 
obtain the photon fiux normalization [22] , the statistical 
uncertainties associated with the photon flux normaliza- 
tion are always far below 1%. For this reason, an overall 
1% uncertainty for the statistical error of the normaliza- 
tion is included as a very conservative estimate of this 
uncertainty. 

Since data collection intervals were taken with two dif- 
ferent incident electron beam energies, it is reasonable 
to compare whether any systematic differences arise be- 



tween the sets of data taken at these two energies. To 
make this comparison, a histogram was generated for the 
ratio R, obtained by taking the differential cross sections 
of the 3.115 GeV data set and dividing those values by 
the cross sections obtained in the 2.445 GeV data set, 
for each kinematic bin. A Gaussian was fit to this his- 
togram, with a resulting centroid of 0.996 ± 0.002. From 
this comparison, we estimate the systematic uncertainty 
introduced by combining data sets from the two electron 
energies to be 0.4%. A fiuctuation in the photon flux for 
differing endpoint energies could be due to the energy 
calibration of the tagger (including effects due to possi- 
ble, slight differences in the tagger magnet field map). 
Since an increase of 27% in initial electron energy (2.445 
GeV to 3.115 GeV) causes only a 0.4% standard devi- 
ation in the photon fiux, a conservative estimate of 1% 
for the systematic uncertainty in photon fiux associated 
with these "field-to-field" differences is used. 

The largest source of uncertainty in the normalization 
is caused by the "tagger efficiency" [6|. The tagger ef- 
ficiency is essentially a measure of the amount of the 
tagged photon beam that survives coUimation, as deter- 
mined during normalization runs. The value of the tagger 
efficiency is dependent upon the electron beam supplied 
by the accelerator, and will vary on a run-by-run basis 
determined by the condition of the electron beam tune. 

To estimate the effects of fiuctuations in the electron 
beam incident on the radiator of the tagger, we calculate 
the proton yield, normalized to incident photon flux, for 
each run, and determine the mean and standard devia- 
tion of this normalized yield (assuming Gaussian statis- 
tics). The run-by-run uncertainty was estimated (see Ta- 
ble [T] ) using the standard deviation of this normalized 
proton yield, dividing by the mean, to get the fractional 
variation. 



E. Absolute normalization 

The systematic uncertainty for the absolute normaliza- 
tion is comprised of flve parts: 

1. uncertainty in the cryogenic target density 
(discussed in detail in Ref. [23|) and length; 

2. statistical error of the normalization (de- 
scribed in the preceding subsection); 

3. run-to-run variations in the normalized pro- 
ton yield (described in the preceding subsec- 
tion); 

4. uncertainty associated with the tagger mag- 
netic field calibration ("field-to-field" differ- 
ences) in combining the results for the two 
incident electron energies (described in the 
preceding subsection); and the 

5. uncertainty of the trigger inefficiency (de- 
scribed in subsection III.H). 
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Table [J shows the vahies for these contributions to the 
systematic uncertainty of the absolute normalization. Af- 
ter adding items (1), (2), and (3) in quadrature and then 
linearly adding the remaining contributions, the system- 
atic uncertainty for the normalization is 3.7% and 3.9% 
for the 2.445 GeV and the 3.115 GeV data sets, respec- 
tively. Even when all of the systematic uncertainties of 
the normalization are simply added together, the result- 
ing uncertainty is 4.4% for the 2.445 GeV data set, and 
4.8% for the 3.115 GeV data set. Since the estimated er- 
ror is never over 5%, the overall systematic uncertainty 
in the absolute normalization is estimated as having a 
conservative upper bound value of 5%. 



TABLE I: Systematic uncertainties in the absolute normal- 
ization. 



Data Set 


Run-to- Run 


Target density 


Target length 


2.445 GeV 


0.98% 


0.14% 


0.28% 


3.115 GeV 


1.25% 


0.28% 


0.28% 


Data Set 


Field-to-field 


Trigger 


Statistical 


2.445 GeV 


< 1% 


< 1% 


< 1% 


3.115 GeV 




< 1% 


< 1% 



V. RESULTS 

The differential cross sections obtained here are 
compared with previous data from MAMI-B ^231, 
GRAAL [H, and CB-ELSA Q in Figs. [7] through [TOl 

For incident photon energies E up to 1.275 GeV, the 
data obtained here are generally in very good agreement 
with previous data. Above that energy discrepancies be- 
tween this work and the results of the CB-ELSA Collabo- 
ration are seen. The disagreement begins at forward me- 
son scattering angles at £^.^=1.325 GeV, and then begins 
to appear at other angles as the energy increases, though 
the differences at larger angles are relatively small. The 
overall systematic uncertainty for the CB-ELSA mea- 
surements is stated to be 5% below 1300 MeV and 15% 
above that energy. This compares with the roughly 5% 
systematic uncertainty obtained here. With these esti- 
mated uncertainties, the data are in statistical agreement 
for the larger angles, while the discrepancies at the small- 
est angles are larger than can be accounted for by sys- 
tematic uncertainties. 

Clearly, additional measurements at forward angles are 
needed to determine whether the rapid increase suggested 
by the CB-ELSA data is correct, or whether the more 
gradual behavior predicted by the fits properly describes 
the trend of the cross section at forward angles. 



VI. AMPLITUDE ANALYSIS OF DATA 

Multipole amplitude analyses provide a powerful tool 
for extracting information about the reaction process in 
as nearly a model- independent manner as possible [1^. 
This approach, in turn, facilitates the identification of 
s-channel resonances involved in the reaction process. 

Cross sections from this experiment have been included 
in, and excluded from, a number of multipole fits to the 
full SAID database [2C)|| in order to gauge their influ- 
ence and compatibility with previous measurements. In 
Table |TT1 the values of resulting from these fits are 
compared to predictions from MAID (MAID05 [2^ and 
MAID03 [13) and an earlier SAID analysis (SM02) 
For the purposes of this discussion, several combinations 
of data sets and fits are presented. A fit called "FDX6" 
was determined from the world database such that it in- 
cludes recent GRAAL [H] and CB-ELSA 3 data but 
does not include the present CLAS dataset, whereas the 
"FD16" fit includes the CLAS data. To emphasize the 
effects of the CLAS data reported here and to minimize 
the influence of CB-ELSA data, the solution "FA06" is 
a fit that the weight factor for this data was artificially 
increased by a factor of 3. 

The fits FD16 and FA06, despite having different 
weights for the CLAS data, are in good agreement with 
each other; this is not surprising, as the older SM02 fit 
also follows the CLAS data. As seen in Table |TT1 the de- 
scription of data by MAID is significantly poorer than by 
any other fit included in this analysis. 

TABLE ILx^ comparison of fits to 3 GeV, fit SM02 (to 
2.0 GeV) J2p], and two recent Mainz fits, MAID05 [H] and 
MAID03 ^{to 1650 MeV [2^]). See text for details. 



Solution 


Range (MeV) 


X^/Data 


FA06 


3000 


55640/25524 


FD16 


3000 


52196/24008 


FDX6 


3000 


49010/23250 


SM02 


2000 


35297/17571 


MAID05 


1650 


141270/21942 


MAID03 


1650 


486266/21942 



The FA06 and MAID fits are compared to data over the 
range of the MAID analysis in Fig. [T] Above this energy 
limit (1650 MeV), we compare FA06 and the older SM02 
solution in Fig. [T] The deviation of MAID03 from data 
and the SAID fits, above approximately 1 GeV, is known 
to the Mainz group. The (preliminary) MAID05 fit is 
much improved. A version to be made publicly available 
through the MAID website will contain further modifica- 
tions [29|]- The MAID05 and SAID fits are significantly 
different in the most forward bump/dip structure and at 
backward angles. 

The forward region continues to differ most in Fig. [51 
where the fits SM02 and FA06 are compared. At in- 
termediate angles, agreement between the CLAS and 
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180 



FIG. 7: The differential cross section for jp — > yr^p below E^. — 2250 MeV. The ang le shown is the pion center-of-mass scattering 
angle. Solid (short dash-dotted) lines correspond to the SAID FA06 (SM02 [2y|) solution. SM02 curves are shown only for 
E.^ between 1650 and 2000 MeV. Dotted (long dash-dotted) lines give the MAID05 [H] (MAID03 ^) predictions. MAID03 
curves are shown only for E^ between 1050 to 1450 MeV. Experimental data are from the current measurement (filled circles), 
MAMI-B ^24| (open triangles), GRAAL 25] (open squares), and CB-ELSA 4| (open circles). The plotted experimental data 
have been selected from energy bins spanning at most 5 MeV. 



CB-ELSA datasets is quite good. Note that the older 
SM02 fit is in perfect agreement with the most forward 
CB-ELSA measurements, though these data were not in- 
cluded in the fit. The model dependence of this forward 



region is further explored in Fig. [8] Note that the FDX6 
result, including CB-ELSA but not CLAS data, is actu- 
ally in worse agreement with the most forward CB-ELSA 
measurements. The FA06 fit (which includes the CB- 
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ELSA points and the data obtained here) yields results 
that also fall far below the most forward point measured 
by the CB-ELSA collaboration. While only suggestive, 
this observation adds some support to the less rapid in- 
crease in the differential cross sections implied by the 
data obtained here. 




180 



60 120 



180 



FIG. 8: The difFerential cross section for 7p tv p at E^, 
= 1775 and 1925 MeV. The angle shown is the pion center- 
of-mass scattering angle. Solid (long dash-dotted) lines corre- 
sponding to the GW SAID FA06 (SM02 [2(1) solution. Dotted 
lines represent FDX6 results. Experimental data are from the 
current measurement (filled circles) and CB-ELSA 0] (open 
circles) . 
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FIG. 9: Fixed angle excitation functions for 7p yr^p. The 
angle shown is the pion center-of-mass scattering angle. Solid 
(long dash-dotted) lines corresponding to the GW SAID FA06 
(SM02 20]) solution. Experimental data are from the current 
measurement (filled circles) and CB-ELSA |3| (open circles). 



The present fits have been generated using the most re- 
cent GW analysis of pion-nucleon scattering data, which 
was extended to cover the full resonance region [30] ■ The 
upper limit on the photon energy is 3 GeV (as compared 
to 2 GeV for SM02 '^]), incorporating the fuU CLAS 
dataset determined here. In Figs. M and [TUl we display 
the energy dependence of cross sections at fixed angle. 
This view most clearly shows the kinematic region and 
scale of disagreement between the CLAS and CB-ELSA 
datasets. Also apparent is the transition, between 2 and 
3 GeV, to a region where the cross sections have a nearly 
energy-independent structure. 

VII. RESONANCE COUPLINGS 

Multipoles from the FA06 fit are compared to the ear- 
lier SAID (SM02) determinations in Figs. [Jl] and [TH 
As suggested by the falling cross sections, the multipoles 



are either flat or rapidly decreasing in magnitude at the 
upper energy limit (e.g., the -©2-^ ^Lud E'^^ multipoles). 
The extension to higher energies provided by this data 
set has resulting in a smoothing of some structures found 
in the 2 GeV limit of the SM02 solution (see, in partic- 

3/2 3/2 

ular, the i?2_ and i?2+ multipoles). Comparisons with 
the MAID05 solution are given in Figs. [Ol and [HI 

Resonance contributions have been fit and the result- 
ing helicity amplitudes are presented in Table IIIII Values 
for the resonance mass Wr, width F, and branching frac- 
tion (F^Tv/F) for the various resonances were taken from 
the recent SAID analysis of pion-nucleon elastic scatter- 
ing [30*]. The electromagnetic resonance couplings were 
extracted using a form 

B{W)il + tT^N)+TBwe"^, (1) 

wherein T^^n was the associated full pion-nucleon T- 
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FIG. 10: Fixed angle excitation functions for yp — > n p below 
1/ib/sr. Notation as in Fig. |9] 

matrix and Tbw was a Breit-Wigner parametrization of 
the resonance contribution. This is similar to the MAID 
form, but allows for a more flexible polynomial function 
to account for the non-resonant contributions from chan- 
nels such as ttA, as described in Ref. [Slj . 

The uncertainties for the helicity amplitudes given in 
Table IIIII extracted in the FA06 analysis correspond to 
a increase of 9 in the fit. This uncertainty is thus 
purely statistical and does not account for subjectivity 
in the resonance extraction and selection of the energy 
range used in the resonance fit. Table IIVI compares 
values from FA06 and the resonance fits over correspond- 
ing energy ranges. 

Most of the resonance couplings determined in this 
analysis generally are in fair agreement with the PDG 
averages, but there are significant disagreements for the 
7V(1650) and iV{1720). The iV(1650) couplings are par- 
ticularly difficult to extract as there is an overlapping 
resonance (the iV(1535) resonance) below this state and 
possibly a third resonance slightly higher in energy. The 
A^(1720) coupling, quoted by the PDG, is not clearly de- 



TABLE III: Resonance parameters for N* and A* from the 
SAID fit to the ttN data [s^l (second column) and helicity 
amplitudes A1/2 and (in [(GeV)"^''^ * 10~^] units) from 
the FA06 solution (first row) and average values from the 
PDG06 m (second row). 
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termined even in sign. Furthermore, the multipoles asso- 
ciated with this state have also changed dramatically in 
the extension to 3 GeV, as can be seen in Fig. [TlJ 

The good agreement between the fit and PDG param- 
eters seen for the A^(1535) resonance also deserves some 
comment. The large PDG error band was given mainly 
to account for the spread in determinations from pion 
and eta photoproduction analyses. This discrepancy has 
largely disappeared [13]. However, the same qualifica- 
tions addressed to the 7V(1650) resonance apply here as 
well. Model dependence in this extraction is certainly 
larger than the statistical error quoted in Table IIIII 
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FIG. 11: Multipole amplitudes from threshold to = 3 GeV for isospin 1/2. Solid (dashed) curves give the real (imaginary) 
parts of amplitudes corresponding to the SAID FA06 solution. The previous SAID SM02 solution [2y] is given by long dash- 
dotted (real part) and short dash-dotted (imaginary part) lines. Vertical arrows indicate Wr and horizontal bars show full V 
and partial widths for F^jv associated with the SAID -kN solution SP06 [30| . 
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TABLE IV: Comparison of the resonance fit and FA06 against 
data for the isospin 1/2 and 3/2 baryon resonances (see text 
and associated Table IIII|I . The "Data" corresponds to the 
number of experimental data points in the given range of the 
Wmin and Wmax limits. Energy limits are those used in the 
most recent nN resonance extractions [3011. 



Resonance 






Fit 


FA06 


Data 




(MeV) 


(MeV) 








iV(1535)5'ii 


1490 


1590 


7129 


7546 


3552 


iV(1650)5'ii 


1620 


1770 


6013 


6520 


2588 


iV(1440)Pii 


1350 


1550 


13434 


13752 


6897 


A^(1720)Pi3 


1650 


1790 


5520 


5717 


2514 


Af(1520)Di3 


1480 


1560 


7252 


7554 


3481 


Ar(1675)Di5 


1610 


1730 


5900 


6139 


2406 


Af(1680)i^i5 


1620 


1730 


5359 


5625 


2187 


A(1620)S'3i 


1570 


1680 


5764 


6018 


2475 


A(1232)P33 


1180 


1270 


6303 


6630 


3351 


A(1700)D33 


1550 


1750 


7285 


7577 


3064 


A(1905)F35 


1770 


1920 


3930 


4089 


1787 



Given the smooth behavior exhibited by the excita- 
tion functions in Figs. [H] and [TOl the CLAS cross sec- 
tions provide no hint of "missing" resonance structure 
between 2 and 3 GeV. The SAID fits implicitly con- 
tain only those resonances found in the corresponding 
SAID analysis of elastic pion-nucleon scattering data. No 
change in the form of the SAID photoproduction fit was 
found to be necessary. In contrast, the CB-ELSA fit 
required many additional resonance contributions, some 
of which are 1- and 2-star rated PDG states, as well as 
a new N(2070) resonance. One possible explanation is 
apparent in Fig. [TOl which shows the CLAS data to be 
somewhat smoother than the CB-ELSA excitation func- 
tions. Model-dependence in the separation of resonance 
and background contributions is also a critical factor. 
This uncertainty can be reduced through measurements 
of further (polarized) data. 



ments. 

Multipole analyses incorporating the CLAS data and 
extending to 3 GeV now appear more smooth near the 
2 GeV limit of previous fits, tending to be relatively fea- 
tureless and rapidly decreasing above this energy. At 
energies above 2 GeV and outside of the angular range 
constrained by CLAS data, the fits were found to be un- 
stable. More progress will require precise measurements 
in that angular region, and also the availability of po- 
larization observables of a similar quality and coverage. 
Corresponding data from nTr"*" photoproduction are also 
clearly required to make an isospin decomposition. A 
partial wave analysis performed in this work using only 
"4-star" resonances at or below 1^=1910 MeV satisfac- 
torily described the data. 
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VIII. CONCLUSION 



The differential cross sections for the reaction 7 p — ^ 
p tt" for incident photon energies between 0.675 and 
2.875 GeV have been measured. These measurements 
have been compared to existing data, mainly from the 
CB-ELSA Collaboration. The overall agreement is excel- 
lent, though statistically significant disagreement is evi- 
dent in some cases between that data and a few values 
obtained here in the forward direction. The variation of 
fits in the forward region suggest further measurements 
at those angles would be very useful. As shown in Figs. [5] 
and [lOl both the SAID fit SM02, developed prior to the 
publication of both the CB-ELSA and CLAS data, and 
FA06 (including both sets) appear to favor the most for- 
ward CLAS measurements over the CB-ELSA measure- 
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FIG. 12: Multipole amplitudes from threshold to i?-y = 3 GeV for isospin 3/2. Notation as in Fig. 1111 
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FIG. 13: Multipole amplitudes from threshold to S.^ = 3 GeV 
parts of amplitudes corresponding to the SAID FA06 solution 
part) and short dash-dotted (imaginary part) lines. 



for isospin 1/2. Solid (dashed) curves give the real (imaginary) 
. The MAID05 solution [2^ is given by long dash-dotted (real 
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FIG. 14: Multipole amplitudes from threshold to E-^ — 3 GeV for isospin 3/2. Notation as in Fig. 1131 
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